clear index
k = 1;
for n = 1:12,
    for m = 0:n,
   
		index(k) = (n * (n + 1) / 2 + m);
        fprintf('index = %d, n = %d, m=%d\n',index(k),n,m);
        k=k+1;
   
    end;
end;

data = load('c:\manoj\projects\WMM2010\GeoidVsNoGeoid\F_Geoid_2009.txt');
FYG = reshape(data(:,5),[360,179]);
figure
axesm eckert4; %initializes projection
framem; gridm;
axis off
set(gcf, 'Renderer', 'ZBuffer'); %%IMPORTANT WITHOUT THIS MATLAB CRASHES
geoshow(FYG',[1,89,0],'DisplayType', 'texturemap')
geoshow('landareas.shp','FaceColor','none');
